function plotBindingPosition(systemData,pivotRotations,hh,pstyle)

if nargin < 3
    hh = figure;
    pstyle = 'b';
elseif nargin <4
    pstyle = 'b';
    figure(hh); hold all;
else
    figure(hh);
    hold all;
end

% [heelPos,bootMoment] = feval(systemData.bootModel,pivotRotations);

forefootratio = 0.346; % fraction of forefoot (foreward of bellows)
% length to total boot sole length

alphaB = pivotRotations(2); % angle from ref for bellows
alphaT = pivotRotations(1); % angle from ref for toe of boot

Lt = systemData.boot.soleLength; %total sole length
Lf = Lt*forefootratio; %forefoot boot sole length (forward of bellows)
Lr = Lt*(1-forefootratio); % length from bellows rearward

initialXheel = [Lt,0]'; %initial heel position
systemData.boot.initialXheel = initialXheel;

xbellows = [Lf*cos(alphaT);Lf*sin(alphaT)];
xheel = xbellows+[Lr*cos(alphaT+alphaB); Lr*sin(alphaT+alphaB)];

xtb = [0,xbellows(1)];
ytb = [0,xbellows(2)];

xbh = [xbellows(1),xheel(1)];
ybh = [xbellows(2),xheel(2)];

Fb = bindingModels(xheel,systemData);

Fb = Fb./norm(Fb);

Fb = Fb.*Lr;

xfb = [xheel(1),xheel(1)+Fb(1)];
yfb = [xheel(2),xheel(2)+Fb(2)];

pstyle0 = 'r';

plot([0,Lt],[0,0],pstyle0,'LineWidth',3);
plot(xtb,ytb,pstyle,xbh,ybh,pstyle,'LineWidth',3);
plot(xfb,yfb,pstyle)

legend('default position','Forefoot','Rear foot','binding force orientation')

